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The Vlasov equation has been solved for the plasma resonance spectra of a realistic model by the 
conductivity kernel method. The results give a clearer picture of the nature of plasma resonance than 
heretofore available. The calculation includes Landau damping and does not, impose unphysical 
boundary conditions. The problem of chief interest here is the scattering of electromagnetic radia- 
tion near the electron plasma frequency from a cylindrical plasma when the wave vector and polari- 
zation are perpendicular to the cylinder axis. The plasma is nonuniform and hounded by a sheath, 
and has diameter a small compared to the free space wavelength. The scattering resonances at the 
lower frequencies are produced by charge density perturbations concentrated at relatively large radii. 
But, the problem of a perfectly collisionless cylindrical plasma cannot be reduced to a one-dimensional 
problem without neglecting some of the resonances with periodic electron orbits. It is argued that 
weak coulomb collisions destroy these "transit time resonances," and that the problem is adequately 
described by keeping just one period of the electron orbit in the calculation of the conductivity kernel. 
The cylindrical problem then reduces to the problem of the steady, driven, oscillations of a thin one- 
dimensional slab of collisionless, Maxwellian, plasma, with a wall at x = which emits electrons and 
absorbs all electrons that return to x = 0. and an insulated wall at x = x w which also absorbs electrons. 
A model is used in which the unperturbed electric field, everywhere in the positives-direction, is uni- 
form in the plasma, 0^X=Ss, and joins smoothly to a harmonic oscillator field in the sheath, s =£ X =S £ u ,. 
The conductivity kernels for a large number of frequencies have been calculated, and inverted, on 
large electronic computers. The results show that much of the Landau damping which determines 
the line shapes is concentrated near the sheath, and the resonance frequencies are determined by 
the properties of the sheath and the neighboring regions of the plasma. 



1. Introduction 

"Plasma resonance" is the name given by Tonks 
[1931a, b] to the resonances he found near the plasma 
frequency in a system consisting of the cylindrical 
positive column of a low-pressure mercury arc placed 
between a pair of capacitor plates which were driven 
by an oscillator, when the plasma axis was perpendi- 
cular to the electric field. Plasma resonance has since 
been found in free-space scattering of microwaves 
from the positive column [Romell, 1951], in backscatter 
of 55 megacycle radiation from meteor trails [Billam 
and Browne, 1955], in the reflection from the positive 
column placed across a waveguide [Dattner, 1957], 
and in noise radiation [Lustig, 1964]. Our purpose in 
this paper is to obtain a reasonable theoretical under- 
standing of plasma resonance in the absence of a 
magnetic field. We calculate the response of a nearly 
collisionless plasma to an externally applied weak 
driving field for frequencies on the order of the mean 
plasma frequency. We include the nonuniformity of 
the plasma and the sheath and do not introduce arbi- 
trary boundary conditions. 

Tonks attempted to explain the strongest resonance 
as follows: Consider a sharply bounded cylinder of 
cold, neutral, uniform plasma in free space. Let N 
be the electron density, m the electron mass, and e 



1 Most of this work was done at the University of California at San Diego, and was in- 
cluded in the author's Ph.D. thesis. 



the magnitude of t he electron charge. A slight uni- 
form displacement, 8x, of the electrons with respect 
to the heavy and therefore nearly stationary ions 
produces a surface charge density — Ne&x cos ^ ar, d 
a uniform restoring force inside the plasma of —277 
Ne 2 8x- The system should oscillate at the angular 
frequency oj, with co = co /J /V2, where co p = (47r/Ve 2 /ra) 1 / 2 
is the plasma frequency, if the internal motion of the 
surface charge is ignored. Tonks identified his 
strongest resonance with the response of this dipole 
plasma mode to the driving field. The observed 
angular frequency was near co p /v2. 

An equivalent model, an infinitely long cylinder, 
which was thin compared to X = 27rc/o), and had dielec- 
tric constant 1 — co 2 /a) 2 , was used to predict the back- 
scatter from meteor trails [Herlofson, 1951]. The 
experimental results for meteor trails [Billam and 
Browne, 1955] are quite interesting. The authors 
were apparently only interested in the main peak, 
but their results seem to indicate that there may be 
as many as six resonances, while the theory predicted 
one. Tonks [1931a, b] suggested that the one or two 
"anomalous" resonances he observed were caused by 
the nonuniform electron density and nonzero plasma 
temperature. Theoretical work since Herlofson's 
paper [1951] has been directed toward taking adequate 
account of these two conditions. The theory [Herlof- 
son, 1951] also predicted that the scattering width at 
resonance would be 8A/77, where the scattering width 
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DISCHARGE CURRENT, I 



FIGURE 1. Microwave reflection from a positive column placed 
across a waveguide; taken from Dattner [1957]. 



was defined as the diameter of a cylinder which would 
absorb and then reradiate isotropically all of the radia- 
tion incident on it, and produce the same backscat- 
tered intensity as the actual scatterer. Romell 
[1951] found that at \ = 30 cm, a 3.2 cm diameter 
cylindrical discharge plasma gave backscatter intensity 
at the first resonance which was 0.8 that of a long flat 
copper strip which was 15 cm wide. At least the first 
two resonances excited by plane waves are essentially 
pure dipole modes [Boley, 1958], as required by the 
theory. Higher multipoles have been excited by 
other driving field configurations [Parker, Nickel, and 
Gould, 1964]. 

Dattner [1957, 1963] has experimented with the 
cylindrical positive column of a mercury arc placed 
across a waveguide with its axis perpendicular to the 
£"-vector. He has demonstrated that there are often 
six or seven resonances (fig. 1). The strongest falls 
at the highest arc current (lowest o>/co p ), and they form 
a regular sequence. Dattner also introduced the 
concept of the "series limit," i.e., a definite upper limit 
on the resonance frequencies. 

The steady state structure of the positive column 
and sheath is only beginning to become well under- 
stood [Ott, 1963; Self, 1963]. Where characteristics 
of the steady arc are needed we shall use experimental 
data [Gabor, Ash, and Dracott. 1955; Gierke, Ott, 
and Schwirzke, 1961; Harp and Kino, 1963]. 

The kinetic theory of plasma resonance is developed, 
in section 2, by the conductivity kernel method. For 
comparison with the kinetic theory, the adiabatic 
fluid equations are solved in section 3. The fluid 
theory is unsatisfactory for calculating the electric 
fields. Nevertheless, with a carefully chosen boundary 
condition, it is capable of predicting the resonance 
frequencies. 

2. Kinetic Theory of Plasma Resonance 

Drummond, Gerwin, and Springer [1961] have de- 
rived a general integrodifferential equation for weak, 
oscillating electric fields in a collisionless plasma. 
We give a simple derivation of the conductivity kernel 
equation, for cases where the perturbed magnetic 
field may be neglected, in section 1. Sections 2-2.5, 

describe the simplifications which lead from the 
cylinder problem to the one-dimensional problem that 



we have solved, and sections 2.6 and 2.7 are descrip- 
tions of the numerical solution and the results. 



2.1. Conductivity Kernel Equation 
Maxwell's equations, 

V B = 0, V-E = 477p, 

cV X B = 4ttj + dE/dt, 



and 



(1) 



cV xE=-dB/dt, 



with the usual symbols, (Gaussian units except /, in 
esu; i.e., j = Nev), imply the general wave equation 



V 2 E-^ 
c 



i a 2 E_47r dj 



c 2 dt 



-7-h47rVp- 



(2) 



Conductors and dielectrics external to the plasma 
could be included by using the constitutive equations 
[Panofsky and Phillips, 1962], but have no fundamental 
role in plasma resonance. 

The source terms are related to the distribution 
function F s of species s (where s = ± stands for singly 
charged ions or electrons respectively) by 



} = e^s \F s (\)d 3 v, 

g J - oc 



(3) 



F s {\)d 3 v. 



The identity 

F f (x, v, f) = F,(xo, vo, k)+ f ^ (x ';, v '»'' 

Jt dt 



dt\ (4) 



where the set of points (x' (x, v, t'), v' (x, v, t')) lie on 
any continuous path S through (x, v), provides the con- 
nection with the kinetic equation, 



dFs 

dt 



dF s 
dt 



(5) 



where S must now be the mean path in phase space of 
the particles in a small neighborhood around (x',v'). 
Consider small, oscillatory perturbations about the 
steady state, 



Fs^NoJos+fse-^, 



^here 



and 



/: 



/os(x, \)d 3 v = 1 , 



(6) 



A=Ao + A,€- ia,t , 
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where /f = E, B, j, or p, and where N 0s {x) is the un- 
perturbed density, and where o> is of the order of the 
electron plasma frequency. 

In the discharges where resonance has been ob- 
served, the frequency of momentum change by colli- 
sions is very small compared to the plasma frequency 
so we tentatively write the linearized Vlasov equation 
for/., 



is zero beyond the turning point, and where we must 
take e~ uoit '~ n to be nonsingular at the off-orbit points 
where it is undefined. Also define 



o- p i = Vi(x,y; x',p)/t>/(x,v; x\ 1), 



(11) 



§=-— E, • V„/o*(x, yW*(x). 
at m s 



(7) 



Collisions acutally play an important role, which is 
described in section 2.7. We have neglected the 
vXBi/c term on the right-hand side of (7), which is 
valid for high frequency modes or for isotropic /o*. 

The zero order part of the set of equations (2) to (5) 
is complicated by the collision terms. We assume 
solutions of this set, and use experimental datafGabor. 
Ash, and Dracott, 1955; Gierke, Ott, and Schwirzke, 
1961; Harp and Kino, 1963] when numbers are needed. 
This process is consistent only if no high frequency 
phenomena play an essential part in determining the 
measured "steady-state" [Gabor, Ash, and Dracott, 
1955; Ott, 1963]. 

Substituting eq (7) into (4), eq (4) into (3) and eq (3) 
into (2), we obtain the following linearized equation <^~(x; x') = V 
for Ei(x): . 8 , p 



which is either of ±1. Now write the orbit integral 
as a space integral 

J Ei(x') • V r'/i),-(v V W ' ©*(X 'W 

= ]To>- [Ei(x')- V,/o,(vV^'ft|(x')^cPjc', (12) 



Vi 



where the integral extends over all regions which are 
reached by any electrons, i.e., the limits of the space 
integration are independent of velocity. 

Using (12) in (8) and exchanging the order of inte- 
gration gives 



V 2 E, + -^-E, 



^-(x; x') • E,(x')dV, (13) 



where 



V 2 E,+^E, 



-""2 



V T V) 

c 



/: 



E, • V i//o*(x\ v')e- iw <' a)'i(x')dt'd 3 v, (8) 



where a>f = 4>7rNo 8 e 2 lms. The quantities x'(x, v, t'), 
v'(x, v, t') are the position and velocity at time t' on 
the orbit that goes through (x, v) at t' = t. Equation (8) 
assumes the initial conditions 



o>(V - l -§ v^-W) ^vfos~ (14) 
c z ~" v. 

Equation (13) is the conductivity kernel equation 
[Drummond, Gerwin, and Springer, 1961 J. 

If the plasma thickness, a, in the direction of the 
dominant fields is very small, i.e., ka<<l, where 
k = co/c, then we may describe the processes inside 
the plasma by the static approximation, V X Ei 
— 0, Bi — 0, most easily obtained by setting l/c = in 
Maxwell's equations. Then, since 



/s(x,v, t - 



-°°) = 0, 



(9) 



have 



so we must follow the Landau convention [Landau, 
1946; Thompson, 1962] with respect to any poles in 
the integrals in (8). Let Xj(x) be a particular direction 
such that v'j # except at isolated turning points, and 
write dt' as dx[\v[. Let /; be the number of times 
(counting backward in time from /) that the orbit 
through (x, v, t) has passed through a particular x\. 
The path integral along each passage of the orbit is 
expressed as a space integral by writing the integral 
as a set of delta functions located on the orbit, i.e., let 



yxvxE, = 0, 



V 2 E! = VV • Ei, 



$ 



iio(t'-t)- 



-i<o(t'-t) 



U{x'j—Xj,(i,)8(x' 



(10) 



where x^ is in the plane perpendicular to x[ at the 

position x' on the orbit, Xj.tp is the turning point, and 
u(x[ — %i f t p) is a step function such that the amplitude 



which, when combined with Poisson's equation, 
indicates that the wave equation, (2), splits into two 
separate equations, 

V 2 E!=47T V pi, 



dEi/d*=— 477ji, 

Equation (13) splits similarly, the more useful equation 
being 

E,(x)= jy (x; x') • £i(x')dV, (13') 
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where 

y(x;x') = 



7?./: 



^-,e7 M ''-' ) w?(x')V,./o s ^. 

(14') 

If Ei = E p + Erf, where Ed is a driving field with sources 
far outside the plasma, only E ; , can satisfy the electro- 
static approximation, and since E r / satisfies the free 
space wave equation, (13') becomes 



f </> .(E p + E d )dV. 



The electrostatic approximation is valid whenever 
the characteristic length is short compared to the 
free-space wavelength. 



2.2. Positive Column, co^co_ 

In the positive column the electron distribution 
function / i, is observed to be very nearly Maxwellian 
[Langmuir. 1925; Gabor, Ash, Dracott, 1955J so that 
fo-O)- is invariant on the unperturbed orbits, which 
simplifies (14). Henceforth we consider only high 
frequencies, 



a) ~ co_ > > au 



(15) 



and neglect the ion contribution, ^~+. Then with the 
Maxwellian distribution 



/„_ = (/3/77pe-^ 2 , 
the conductivity kernel becomes 
^"(x; x')=-2a> 2 (x)2: J <Tpi 

„ J -00 



(16) 



▼ — T v ) ^ Mr '-°/3/o-(v) >«, (17) 



and in the electrostatic limit the kernel is 
.x (x; x ) — > dpi — e i xm "fo-ivjcPv. 



(17' 



2.3. General Method of Numerical Solution 

In the above form, (13) corresponds to a boundary 
value problem. The problem can be changed to a 
system of linear equations, which in principle can 
always be solved numerically. In the experiments 
the externally supplied driving field, E^, is always 
known. Split E into 



Ei — E r / + E ;> , 



(18) 



V 2 E„ + ^E„ = 



^(x;x')(E ; ,(x') + £ rt (x') )dH', 



(19) 



since E^ satisfies the free-space wave equation. 

In matrix notation (19) has the form A u = b, with 
solution y—A~ x b, where 

y= column (y u y 2 , ys), 

yj = column (£ PJ <Xi), £ P j(x 2 ), . . .), 

where xi, x 2 , . . . exhaust the space of x. The in- 
homogeneous terms, 6, is constructed from the driv- 
ing term in (19) and has the same labeling arrangement 
as y. The matrix A is a 3 X 3 supermatrix 

A = [A U ], 

where the matrices Aij are, with 1, m, being the labels 
given to points Xi, x m , and i, j being the labels on the 
components of E, 



[A,] 



-8ij 



~ 7 6iw"t"^>l; 

c z 



\-ffij\lm, 



where E,, is due to sources in the plasma. Then (13) 



where LimEj is a difference representation of V 2 £j. 

In principle A and A 1 can be formed numerically 
for any system, e.g., cylindrical geometry. But in 
the cylinder problem (with no z-dependence), c^has 
four components, each of which is a complex array 
formed by integration over two-dimensional orbits 
and is expensive in computer time. Inverting A is 
also expensive, since it would necessarily be quite 
large. In the following two subsections we reduce the 
problem of the low-lying dipole modes of a thin cylin- 
drical plasma to a one-dimensional problem, which 
we have solved numerically. 

2.4. Transit Resonance and the Coherence Length 
of an Orbit 

Drummond, Gerwin, and Springer [1961] pointed 
out that the orbit integral could have poles which arose 
from the resonances between the orbits and the fields. 
These transit poles arise whenever oj/27t is a harmonic 
of the transit frequency 1/7, where T is the period of 
the particle orbit. For a slab of width L, with zero 
magnetic field and reflecting walls, T=2L/v. In a 
magnetic field there are transit time poles at the cyclo- 
tron harmonics, and their effects have appeared in 
experiments [Bekefi, Coccoli, Hooper, and Buchs- 
baum, 1962]. However, they do not appear in the 
ionosphere backscatter experiments [Bowles, private 
communication]. The apparent reason has been given 
by Dougherty [1964], who solved a model Fokker- 
Planck equation, and found that the leading effect 
is not a simple collisional damping [Comisar, 1963], 
but rather a term which tends to destroy the poles. 
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Consider a collisionless system in which all orbits 
are periodic. Go back to (12) and sum over successive 
periods of the orbit, which are identical except for a 
phase factor e ipajTo , where 7o(x, v) is the transit time on 
the collisionless orbit. Then the orbit integral becomes 



J' Ei(x') _ Vv>Mv>- i(ot '(o 2 s (x')dt' 



Ei(x') • Vvfose-^ix'W, (20) 



which imply that 



= U> in 



where the integral includes one complete period of 
the orbit. The transit time poles arise from the sum 



£ e ina>r 0=1 /( 1 _ e io>T ) > 



(21) 



Since we have already adopted the Landau [1946] 
convention, we may assume an infinitesimal positive 
imaginary part of a> to force this sum to converge to 
(21). The transit time poles are at 



2irm\T(u m = 0, 1, 2, . 



(22) 



Note that the conductivity kernel equation becomes 
exact if we integrate along exact unperturbed orbits 
including collisions. The conductivity kernel equa- 
tion then becomes equivalent to the single particle 
Liouville equation. Each orbit integral may be re- 
placed by an average over many trials of the orbit. 
We assume it is sufficient to include only random, 
weak, coulomb collisions, so that the exact orbits 
are only slightly different from the collisionless orbits 
over times short compared to the relaxation time, l/i^o. 
Then the magnitude of the left-hand side of (20) is 
less than 



2 < e 1 ™ 7 ™ > Wx, v), 



(23) 



where / max ( x > v ) is the set of maximum values attained 
by the orbit integral over one period of the orbit, in- 
cluding points x not reached on the collisionless orbit, 
and Tnc is the transit time on the nth passage of a trial 
of the collisional orbit. We may reasonably expect 
that only the first term in the sum in (23) will be im- 
portant if the TVs have a spread, 7\ T, satisfying 



AT^2tt/oj. 



(24) 



For v on the order of y t h, and VqT < < 1, we have the 
approximate relations 

T v ' 



and 



AT^v r 



for weak collisions. Let T= T c when the equality in 
(24) is satisfied. Then 




If A T is the width of a Gaussian distribution of orbit 



V^- vqTv, 



times, we find that ^ l e inuTnc\ = y e -yn*+ma>T 0t) wne re 

y = ((OvT^2f. When T > T c all but the first term in 
the sum are negligible, which supports (24). 

For a plasma which is thick compared to a Debye 
length, Ld, A T may be much less than T and still 
satisfy (24) for all electrons except those with such 
high energies that their number can be neglected. 
The condition for the existence of transit time poles 
is that the "coherence length" of the orbit must in- 
clude many passages. This requirement is not sat- 
isfied by low-pressure mercury arcs with no or weak 
magnetic field or by the ionosphere [Dougherty, 
1964]. With L c = (n<o P lv ) ll2 L D we find L r =WL D for 
Dattner's data [1963], using standard estimates of the 
momentum changing frequency due to weak coulomb 
collisions, v {) . 

However, we need not abandon the conductivity 
kernel formulation for the Fokker-Planck equation. 
The range R, of the kernel is short compared to typical 
orbit coherence lengths for almost any interesting 
plasma. To zeroth order in RIL ( we may neglect 
collisions during the first passage from x to x' . Since 
the size of the plasma is many Debye lengths, e ,0)Ti) is 
a rapidly oscillating function of v; the contribution to 
the kernel of any given passage after the first is very 
small. Even if L ( is several times the plasma size 
the kernel will not be significantly affected by any 
but the first passage. We simply drop the higher 
passages on the orbit. Thus, while collisions nowhere 
appear explicitly in our conductivity kernel, they play 
a very important role in determining its form, and this 
is in spite of the fact that by the usual criterion, 
vq < < to, collisions may be neglected. 



2.5. Reduction of the Cylinder Problem 

The highest noticeable resonance is generally ob- 
served to fall near the plasma frequency at the center 
of the discharge tube [Dattner, 1963; Schmitt, 1964], 
and the so-called "main" resonance, at roughly 
[Dattner, 1963] 0.4 of the highest plasma frequency 
(we have assumed that the average electron density 
is approximately one-half of the maximum). The 
lowest modes should be concentrated at relatively 
large radii, where co ~ a),,(r). Their character is essen- 
tially one-dimensional for low-multipoles and long 
wavelengths in the z-direction, provided that transit 
time resonances are destroyed by fluctuations. 
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Consider the scattering of plane waves by a thin, 
infinitely long, cylindrical plasma with z-symmetry. 
The driving field is 



E„ = 



O i(ky-o)t) 



or 



Ed=(r cos 6+ § sin 6) [cos (kr sin Q)+i sin (kr sin 6)]e~ i(ot , 

(25) 

where r, #, are cylindrical polar coordinates and x, y 
are rectangular coordinates in the r, 6 plane, with z 
along the cylinder axis and = on the #-axis. 

For a thin plasma of radius a < < A, where A is 
the freespace wavelength, we have kr^ka<<\. 
In Dattner's experiments [1963] ka ranges from 0.08 
to 0.33. Write 



E d =xe- i<at + l 



%a, 



(25') 



and neglect terms of order kr for r ^ a, while retaining 
the full wave equation for r> a. This has been called 
the quasi-static approximation [Kaiser and Closs, 
1952]. To zeroth order in ka the conductivity kernel 
equation is (13'), where the kernel is given by (17') 
without the sum over higher passages: 



~ gjfogfr) f- w 

J^(x; x ) = — r/o-(v)e. 



-M*'-0 3 



(17") 

From the symmetry of the unperturbed plasma the 
^-dependence of S^ at fixed r must be of the form 
g(0' — 0), which is sufficient to guarantee that driven 
oscillations have the same multipolarity as the driving 
field. Therefore, to zeroth order in ka the fields are 
pure dipole modes, 



:<f(r)r cos + e(r)0 sin0, 



(26) 



$f'+ — + -) cos = -4,<ireNi(r) cos 0. 



V E D 



To zeroth order in ka [see the discussion above (13')] 
VXE p = — (€+re'-#)z = 0. 



Let A. r = e/€' be the characteristic length and assume 
(as we have done everywhere in the static approxi- 
mation, and as is verified by the results) that X r /a < <1. 
Then to lowest order in ka and \ r ja 



>, 



,7a). 



Without transit time poles the terms of order \ r ja 
may be neglected when we substitute (26) into (13'), 
i.e., the off-diagonal kernel S r0 does not contribute 
to the equation for $ because e is small. 



Since the range, /?, of the kernel in any direction 
should be only a few Debye lengths, the 6' dependence 
may be expanded, 



cos 6' — cos 0- 



(Rldf. 



Neglecting the term 0(R/af, the cos dependence 
can be cancelled out of the conductivity kernel equa- 
tion for all 6 such that cos 6 > > ka. 

By choosing xi = r we can now do the 6' and z 
integrations immediately [see the definition of #:**#'-') 
(10)], and obtain the one-dimensional equation 



'(&■ 



f V(r; r' 

Jo 



){£ , {r') + E d {r'))dr 



2iBcx) 2 Sr) f x 

X(r; r') = — i/ -(^ M ''-°^, 

co J-.C 



(27) 



(28) 



where e~ laj(r ~^ is just a step function at the turning 
point multiplied by e^^'^K The orbit times, t' — t 9 
are now calculated in a one-dimensional potential, 
the corrections being \R\d). 

For a given energy there are two terms in (28), 
one connecting r and r directly, the other connecting 
them through a reflection from the nearest turning 
point. The kernel is symmetric. 

For r ^ a we have 



V^p + ^Ep^O. 



(29) 



The appropriate solution of (29) must be chosen by 
requiring continuity of E p at r= a. Therefore a meas- 
ure of the scattered radiation is 



e(a) = I dr + e(r ) 



,a — r 



a — r 



(30) 



If the radiation does not penetrate significantly into 
the plasma beyond r , then 



e(a) = \ P g(r)dr, 

V Jr 



(31) 



e(r ) 



so that to zeroth order in — — the radiation fields are 

e(a) 

proportional to the "dipole moment," 



Jr 



r)dr. 



(32) 



The outgoing wave solution of (29) for the ^-component 
of E p is . 



e(r) = AH[ 1 \kr), 



(33) 



where k = co/c and H\ l) is the Hankel function of type 1. 
The boundary condition is 



e(a) = D/a. 



(34) 
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The average power radiated away per unit length is 



Pre 



lim (/€ 2 (/))/4. 



Define the incident power as 

P lM =2a£>/4ir, 

and define the dipole scattering width, d, by 
d = ar ra Jr inc- 



(35) 



(36) 



(37) 



The dipole scattering width is one half of Herlofson's 
isotropic scattering width. For small ka the first 
term in the asymptotic series for H\ l) gives A, and the 
scattering width is given by 



d 
a 



17 (l » DD * 



(38) 



The approximations we have made to reduce the cylin- 
der problem to a one-dimensional problem correspond 
to the following model, in which it is assumed that the 
y, z dependence of Ea can be neglected: The slab of 
plasma is bounded by a flat inner wall at ^ = and a 
flat insulating wall at x — w. The inner wall emits a 
half-Maxwellian distribution of electrons and absorbs 
all electrons which strike it. A sheath is built up on 
the insulating wall. The infinitely massive ions are 
distributed in such a way as to satisfy Poisson's equa- 
tion. This model could be duplicated very closely 
if the inner wall were a hot tungsten plate emitting 
cesium plasma. 

It is convenient to scale (27) and (28) by measuring 
distances in Debye length times vTand time in plasma 
periods. Let x = be an arbitrary position, and let 



L D = (2kTlma$(0)yl 2 . 



The scaling is 



x- 



xLiu 



t — t'-*TlG)p(0), 

<£<)—> — kTipole, 
E-*-kTE/L D e, 

CO -*(i)p(0)ft, 



(39) 



(3v 2 



w, 



where </? (0) = 0. We will usually pick i = at the 
inner wall. The scaled longitudinal kernel becomes 



K(x 9 x) = 1 7= e~ w V e! mx > w > ±5 x>) dw, 



(40) 



where ± is the sign of v. Note that , since <p (x) + w = w 
is invariant on the orbit, the kernel is symmetric. 



The insulating wall which confines the positive 
column is not included in the conductivity kernel as 
written above. Unperturbed orbits with energy above 
the floating potential do not return to the plasma, which 
is accounted for by not including orbits after they 
have hit the wall. The oscillating electric fields will 
cause some electrons to reach the wall even though 
their unperturbed orbits do not, and vice versa. This 
causes a perturbation in the number of reflected orbits 
and therefore in the current and the electric field. 
Let v c (x) be the velocity of an electron whose unper- 
turbed orbit reaches the wall with v' = Q. i.e., 



v' c {x, v c \ x w ) = 0. 



(41) 



To first order in the perturbation 8v c of v c , the per- 
turbed current at (x, t) due to perturbations in the wall 
current is 



8/ = - epN {) (x) ( - 1 e'K U) vAx)8vr(x, t). (42) 
Integrating along the unperturbed orbit 

p ft-t'ix, v c , jr w ) 

8v c = Eiixy-t'-'W 

m J- x 



i: 



_ e i<d>( X , v c , xj E x (x')e itoT dT n (43) 
m 



where the orbit starts at the wall. Because we drop 
all transit time resonances 8v c is nonsingular. There- 
fore 



tj=0er*t 9 



(44) 



which may be neglected if we consider only w's such 
that the oscillations are concentrated in regions where 
exp {— /3iif.{x)) < < I. In the plasma exp ( — /Jl£) 
= 0(mlm + yl 2 . 

A similar one-dimensional problem — driven oscil- 
lations of the sheath on a semi-infinite uniform 
plasma — has been solved numerically by Pavkovich 
[1963] for frequencies co < co t) (s). The properties of 
this system and the non-uniform plasma with a sheath 
are very different. 



2.6. Numerical Solution of the Conductivity Kernel 
Equation 

The numerical solution of the conductivity kernel 
equation was performed on the computers of the 
University of California at San Diego and the Na- 
tional Center for Atmospheric Research. One pro- 
gram formed the conductivity kernels and stored them 
on magnetic tapes. A second program read the ker- 
nels from the tapes, called a standard subroutine which 
inverted the kernels by the Gauss-Jordan method, 
and formed the solution, E p . 

Efficiency of the inversion program required that 
the matrices needed for the inversion be stored en- 
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tirely in memory, which limited us to about 70 points 
on the x-axis. Actually, 61 points were used, and this 
was very satisfactory. For several reasons efficiency 
decreased especially rapidly for errors in the kernel 
of less than a few tenths of one percent, so errors were 
held to that level. The resulting errors in E p range 
from less than one percent to a few percent. The 
calculation required one hour of machine time for a 
batch of 20 values of il. 

All integrations were done as simple sums. When 
the path length, x — x\ was long, the phase, Qt{x, v; x') 
sometimes varied by even more than 2tt between 
neighboring values of v on the integration mesh, so 
that even though the mechanics of the integration were 
unchanged (simple sum) those portions of the velocity 
integration were actually done by sampling [Meyer, 
1954]. Accuracy of the integration was checked both 
by examining the convergence of the kernel and the 
fields as the step sizes were made smaller, and by 
comparison with the Landau kernel [Pavkovich and 
Kino, 1963; Pavkovich, 1963] in the case of very weak 
fields. 



2.7. Results From the Conductivity Kernel Method 

Experimental data [Gabor, Ash, and Dracott, 1955; 
Gierke, Ott, and Schwirzke, 1961; Harp and Kino, 
1963] on arc structure was used as a guide, but be- 
cause of considerable uncertainty in the plasma 
frequency at the sheath edge it was necessary to try 
various values of the steady fields. Since the field 
measurements and plasma resonance experiments 
have not been carried out on the same plasmas, 
there was no point in searching for steady fields 
which would give D(il) closely resembling, say, Dat- 
tner's curves. 

All of the systems studied had a constant field, Eo. 
in the plasma, smoothly joined to a harmonic oscil- 
lator field, (x — Xo)flsh , in the sheath, where xq < s is 
such that Eo = (s — Xo)ilg h . We give results for two 
systems which differ in Eq, ft S h<, and thickness. With 
Eo and ft*/, scaled to the sheath edge, system 1 has 
Eo = 0.121, fish = 0.746, which is a reasonable fit to 
Gabor, Ash. and Dracott [1955] data. System 2 
has stronger fields, £0 = 0.232 and £1^ = 1.93. The 
relative sheath strength, df^/iso, is about three times 
greater than in system 1. System 1 is relatively 
thicker, i.e., over the full range of fl the total oscil- 
lating field has insignificant penetration to the inner 
wall; but system 2 is damped by the inner wall if 
a > 0.40. 

The most important characteristic of the oscillations 
excited in the plasma is the magnitude of the dipole 
moment, 



D m = {DD*yi 2 



(45) 



The energy stored in the macroscopic variables at any 
time is equal to the amplitude of the field energy, 



'-hi?*™*' 



(46) 



The rate of dissipation of energy must balance the 
power absorbed from the external field; i.e., 

p diss = - g^; J () E d ImE„dx, (47) 



or 



P diss =--0£ d ImZ)/&r, 



(48) 



for uniform E ( \. The number of cycles required to 
dissipate the macroscopic energy, S, would be 



Q=$l(2iTE d ImD), 



(49) 



and the phase shift of the scattered radiation is0=tan -1 
(ImD/ReD). Figures 2a and 2b are plots of $ , D m , 
ImZ), Q, versus H for systems 1 and 2, respectively. 
The frequencies are scaled to the plasma frequency at 
the inner wall. Both systems show peaks in the 
dipole moment. 

System 1 is the more interesting. The dipole mo- 
ment and energy have two strong resonances and at 
least tw r o more weak peaks. The strong resonances 
occur at nearly the same frequencies for both <^and D. 

Calculating the scattering width at the first reso- 
nance from (38) we find 

where we have set a equal to the thickness of our slab 
of plasma. At ka = 0.3, which is certainly near the 
upper limit for the quasi-static approximation, d\a= 15, 
or dj\ = 0.1. For &a = 0.1, d/a = 5 and d/\ = 0.08, 
and even A:a = 0.05 gives d/a = 2.7, d/k = 0.02. A 
cylindrical plasma corresponding to system 1 would 
certainly have d\a > 1, and might come within a factor 
of two of Herlofson's cold plasma result, g?/X = 4/7T, 
which was obtained as the radiation damping limit. 

The relative damping, — , is uncorrelated with S, 

but the coupling to the driving field, D, has sharp peaks 
at the resonances. In this sense plasma resonance 
is due to maxima in the coupling rather than the other 
possibility, minima in the damping [Leavens, 1963]. 
Figure 3 shows E p for several values of (1. Beyond 
the second node of E p the extrema are very weak, 
because of Landau damping. The damping may even 
be large for the first half-cycle of E p if it falls near the 
sheath edge, because the electrons reflected from the 
sheath tend to be out of phase with the unreflected 
electrons [H sh ~ 0.75ft p (s)]. This is apparent in the 
kernels, three rows of one of which we show in figure 
4. The two rows with the diagonal far from 5 are simi- 
lar except for a scale factor. 
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FIGURE 2. Macroscopic energy, dipole moment, imaginary part of 
the dipole moment, quality, and phase shift versus 12 for: (a) 
System I; (6) System 2. Scaled to plasma frequency and Debye 
length at X = 0. 

Driving field= I. 



System 2 has one strong resonance and two weak 
resonances. Above (1 = 0.40 the field penetrates to 
the inner wall, which apparently causes enough damp- 
ing to prevent any more peaks from occurring. Multi- 
plying the frequencies in system 2 by two will scale 
them to the same base as those in system 1. The 
positions and spacings of the peaks are quite different 
from system 1. The Q(Cl) curve has even less relation- 
ship to the D(fi) curve. The first resonance is weaker 
than in system 1, but could still have dja > 1 for reason- 
able values of a and k. 

We have studied several other systems with stronger 
steady fields. One of them was a very thin system, 
so that the total E=E p -\- Ed penetrated to the inner wall. 
All showed well defined peaks in />(H), but none had 
the strong resonances of system 1. The effect of a 
stronger field on the electron orbits is to decrease 
t(x, v; x') 9 lowering the minimum velocity for which 
(ot < < 2tt, and thereby enhancing the imaginary 
part of K near the diagonal and weakening the real 
part, so more damping is to be expected with stronger 
fields. 

We would like to emphasize the following: 



Any relationship between the first resonance fre- 
quency and the average plasma frequency is purely 
an accidental consequence of the form of the un- 
perturbed plasma. This is clear because the portions 
of the plasma into which the fields do not penetrate 
can be modified freely without affecting the calcula- 
tion. This implies that the partial success of cold 
plasma theories [Tonks, 1931a, b; Herlofson, 1951: 
Crawford, 1963] in describing the first resonance is 
also accidental. There is no difference of kind 
between the first and higher resonances. It is 
especially worth emphasizing that all of the resonances 
come from a single model, and one should not add an 
extra cold plasma resonance to the resonances found 
from kinetic theory or fluid theory [Crawford, 1963]. 
There is nothing fundamental about the series limit 
[Dattner, 1963; Crawford, 1963J. The higher reson- 
ances produce relatively weak changes in the dipole 
moment because of Landau damping. Furthermore, 
the Landau damping can be expected to increase when 
O is greater than the largest plasma frequency. These 
two effects presumably bring resonances above the 
"series limit" down to an undetectably small ampli- 
tude. 
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3. Adiabatic Fluid Theory 




3 




I ImK 


2 

1 











-1 


\V 1 DIAG0NAL AT IX= 15 


-2 


\ /\/ReK 

i i i i i 



10 



20 



30 
Ix 



40 



50 



60 



FIGURE 4. Three rows of the conductivity kernel. 

ft =.255. System 1. 



m^l + V{N<v2>) + NeE = ^ (51) 



where TV is the electron density, and 

v n >=\ v n f(v)dv. (52) 

J — X 



< 



In this section we, (1) sketch the derivation of a 
linearized fluid-type equation of motion by the usual ln what follows we use v=<v>-fu, < u > = 0, 
procedure of taking moments of the electron Vlasov 



equation and terminating the set by making the 
adiabatic approximation; (2) consider possible "bound- 
ary conditions" on the fluid equation; and (3) calculate 
resonance spectra. 

3.1. Fluid Equations for Electrostatic Oscillations 
The first two moments of the Vlasov equation are 



< v% > = < vl > + < u\ > 



:3<i4>. (53) 



^+V.(7V<v>) 



(50) 



The set of moment equations is terminated by taking 
the second moment equation and setting < muju k > 
= 0. The same approximation is made more easily 

by the adiabatic approximation ( h < v > • V 

(PN~ y ) = 0, with y— 3 since only one degree of free- 
dom is involved. Linearizing the moment equations 
and using Poisson's equation we find the equation of 
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motion for one-dimensional problems, 

2eE' 



3eE 1 



mur 



ir 



-o>%{x)- 



Ep — —-Ed- 
it 1 



The normal form of (54) is 

y"+(q 2 -p 2 )y-- 

where 



o) 



E .E, t 



It may be, however, that we can find approximate 
conditions for resonance which would be used in the 
same way as boundary conditions. The conductivity 

(54) kernel solutions suggest that we might be able to 
predict resonance frequencies with the fluid equations 
by using the E\{s) = boundary condition in the 
problem of free oscillations, but that only the odd eigen- 
values are resonances. For slowly varying amplitude 

(54') (as a function of f!) this corresponds roughly to maxi- 
mum dipole moment. The resonance splitting is 
twice that given by a rigid wall boundary condition 
at the sheath edge. 



< ? 2 =w 2 /« 2 ,p 2 =K + 2 



eE' 



9 (eE \* 
WW)**- 



and 



E p = (o p {x)y. 



(55) 



To derive (54) from the moment equations, terms in 
Vq, V ?;<), and V ir are neglected. For the electrons this 
is justified to within a few kT (potential drop) of the 
insulating wall. The fluid equation scales with the 
Debye length. The same scaling used in section 2 
gives 



E'^+E^+^^-^ + E^E,- 



=lw p E„, 



(56) 



or 



y"+3tn°- 



-0!+— - 



H 



El]y = -il p E d . (56') 



Different fluid equations are sometimes obtained by 
the isothermal approximation or by making the adia- 
batic approximation in the laboratory reference frame. 
The differences can be important. For example, if 
<po is a harmonic oscillator potential the requirement 
that E p be zero at infinity gives a continuous spectrum 
of free {Ed = 0) oscillations for (54), but a discrete 
spectrum for the similar equation obtained by making 
the adiabatic approximation in the laboratory frame 
[Weibel, I960]. The latter result disagrees with 
results obtained from the Vlasov equation [Weibel, 
I960]. It is well known that the adiabatic fluid 
equations can give reasonable values for the splitting 
of the resonance frequencies [Gould, 1960; Leavens 
and Abramoff, 1961; Fejer, 1962; Crawford, 1963; 
Leavens, 1963; Nickel, Parker, and Gould, 1963; 
Weissglas, 1963; Hoh, 1964]. The conductivity kernel 
solutions suggest a "resonance condition" which we 
apply to the fluid equations in the next section. 



3.2. Resonance Condition 

There is no reason to impose boundary conditions 
on the plasma oscillations, and indeed the kinetic 
theory results indicated a continuous spectrum. 



3.3. Resonance Spectrum 

In the domain ^ x ^ s, the coefficient p 2 in (56') 
is relatively slowly varying and for weak fields the Eq 
and El) terms may be neglected, only the variation of 
£lp being significant. The free oscillation equation is 



where 



y+|(n 2 -njw)y=o, 



tjp llpy. 



(57) 



The asymptotic solution of (57) which hecomes 
exponentially small into the plasma is [Heading, 1962| 



r _ il„(x) 



J>) % f(il> 



2 :: 



(58) 



g i (x) = -ai--nf,( X )), g>{xo)=o, 

and Ai{^) is the Airy function. The only zeroes of Ep 
are in Ai(^). Let Ai(—£i) = Q, then the condition for 
resonance is 



3 f x 
2}. r / x)dx 



2/3 



:£, i=l, 3, 5, . . . . (59) 



For constant E this reduces to 

^= 12 = - nivi - tanh-VO; v\ = ( 1 - e^/flf J 1 / 2 . (60) 

Solving (60) graphically for system 1, we find the follow- 
ing frequencies for the first five resonances: 



i=l 3 5 7 9 

6 = 2.34 5.52 7.94 10.04 11.94 



ft = 0.178 0.239 



0.286 



0.324 



0.365 



The f/ were taken from the British Association tables 
of the Airy function [Miller, 1946]. The H, are scaled 
to the plasma frequency at x = 0. The agreement 
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between the H? and the positions of the peaks on the 
D(fl) curve for system 1 is quite good; fli is in error by 
1 percent, and the average spacing is extremely close. 
For system 2, we find 



n,- = 0.103 0.156 



0.194 



0.228 0.259. 



The first resonance agrees extremely well with the 
conductivity kernel solution, but the splitting is a factor 
of two too small. 

It is interesting to compare the fluid solutions for E p 
(fig. 5) and the kinetic solutions, figure 3. The plotted 
fluid solutions were actually obtained by numerical 
integration of (56), but they agree with (58). Direct 
comparison is impossible, one set of solutions being 
for the driven problem, the other for free oscillations. 
Nevertheless, it is clear that the fluid solutions are not 
very satisfactory. The fluid equations are useless 
well into the sheath, where their solutions exponentiate 
extremely rapidly [Leavens, 1963]. Otherwise, the 
main qualitative difference is the more rapid decrease 
of the amplitude of successive extrema as one moves 
toward the wall in the case of the conductivity kernel 
solutions. The success in calculating the resonance 
frequencies merely reflects the fact that the spatial 
nodes of the fluid solutions have approximately the 
correct spacing. 



4. Summary and Conclusions 

Plasma resonance in nonuniform plasmas bounded 
by sheaths has been studied by the conductivity kernel 
method. Although this is a collisionless method we 
were able to account for the dominant effect of weak 
collisions by cutting off orbit integrations when the 
path length became greater than the coherence length, 
i.e., when an ensemble of electrons released from the 
point (x 9 v) would reach x after times which would 
have a spread of more than 27r/co. The coherence 
length is much shorter than the mean free path; for 
Dattner's arcs it is approximately one plasma diameter. 
With the orbit cutoff it was possible to reduce the 
problem of low lying dipole resonances of a cylinder 
to a one-dimensional problem. 

The one-dimensional problem of steady state driven 
oscillations was solved numerically for a model which 
used a linear potential in the plasma, continuously 
joined to a harmonic oscillator potential in the sheath. 
The model has plasma resonances. The most inter- 
esting static potential (which was also probably the 
best fit to experimental measurements of the static 
fields near the sheath) had two strong resonances and 
at least two weaker ones. 

The program to solve the conductivity kernel equa- 
tion will work for any static potential, but as yet reson- 
ance experiments and static field measurements in 
the sheath and the plasma have not been done on the 
same arc. Cesium plasmas would be ideal for such 
an experiment. With a hot, flat, tungsten plate facing 
an insulating plate across a gap of ~ 50 Debye lengths 




FIGURE 5. Fluid solutions of the free oscillation problem, E p 
versus position. 

System 1. 



the experiment would correspond exactly to the one- 
dimensional problem solved by our program, and the 
field measurements would not be too difficult. 

Although Landau damping of the conductivity kernel 
solutions was important, we found that the adiabatic 
fluid equations, with a ''boundary condition" chosen 
to correspond approximately to maximum coupling to 
the driving field, gave an extremely accurate prediction 
of the resonance frequencies in the most realistic 
case studied, and predicted the first resonance quite 
accurately (i.e., agreeing with the conductivity kernel 
solutions) even for a system with overly strong static 
fields. Qualitative comparison with Dattner's data 
[1963] shows that the relative spacing of our calculated 
resonances is too large for the higher resonances, 
which is expected because we have used a model which 
has an increasing density gradient into the plasma. 
The "level spacing" is, however, in qualitative agree- 
ment with Dattner's data. 

In conclusion, we note that the "boundary condition" 
was inferred from the conductivity kernel solutions in 
the neighborhood of the sheath. The sheath not only 
influences the line shapes through Landau damping, 
but also has a strong influence on the positions of 
the resonances. 
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